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We develop a semiclassical theory of Coulomb blockade peak heights in chaotic quantum 
dots. Using Berry's conjecture, we calculate the peak height distributions and the correlation 
functions. We demonstrate that the corrections to the corresponding results of the standard 
statistical theory are non-universal and can be expressed in terms of the classical periodic orbits 
of the dot that are well coupled to the leads. The main effect is an oscillatory dependence 
of the peak heights on any parameter which is varied; it is substantial for both symmetric 
and asymmetric lead placement. Surprisingly, these dynamical effects do not influence the full 
distribution of peak heights, but are clearly seen in the correlation function or power spectrum. 
For non-zero temperature, the correlation function obtained theoretically is in good agreement 
with that measured experimentally. 

PACS numbers: 73.23.Hk, 05.45.Mt, 73.20.Dx, 73.40. Gk 

I. INTRODUCTION 

The Coulomb blockade is a fundamentally classical effect in microstructures — the addition of an electron to 
an isolated microstructure requires a certain amount of electrostatic energy, the charging energy e 2 /2C where 
C is the capacitance of the structure. It is the simplest effect of electron charge in microstructures and has 
been extensively studied with regard to both fundamentals and applications in single electron transistors Qj. 
One common way to study the Coulomb blockade is by measuring the conductance through a nearly isolated 
nanoparticle (using tunneling contacts) as a function of a gate voltage which tunes the electrostatic potential 
of the particle. For most values of the gate voltage, the conductance is very small since the flow of electrons is 
blocked because the charging energy is not available. However, when the gate voltage is tuned so that states 
differing by one charge have the same energy, there is a peak in the conductance. The height of this peak 
is simply the conductance of the two tunnel barriers in series, and the spacing of the peaks is uniform with 
separation e 2 /C. 

For the smallest quantum dots and at low temperature, however, quantum mechanical interference becomes 
important. Interference causes variation in both the height and spacing of the conductance peaks. For the 
spacing, single particle quantization and the residual interactions among the electrons are important. For the 
height, the nature of the wave functions become critical: if the wave function of the state at the chemical 
potential is poorly coupled to the leads — if it has nodes at the leads — then the conductance peak is small, but if 
the wave function is well coupled to the leads then the peak is large. In this paper, we restrict our attention to 
fluctuations in the conductance peak heights and investigate what this tells us about wave functions in quantum 
dots. 

Since dots are generally irregular in shape, the classical dynamics of the electrons is chaotic, and so the 
characteristics of Coulomb blockade peaks reflect those of wave functions in chaotic systems [|]-f|] ■ Previously, a 
statistical theory for the peaks was developed ^|J|] by assuming these wave functions to be completely random 
and uncorrelated with each other. The random matrix theory used was known to be a good description of 
energy level statistics, and so likely to be reasonable for wave functions. The experimental data ^,|| for the 
distribution of the Coulomb blockade peak heights were found to be in excellent agreement with the predictions 
of the statistical theory, thus supporting the conjecture of effective "randomness" of the quantum dot wave 
functions. 

A potential problem with the statistical theory was, however, evident in one of the first experiments: there is 
no correlation between different wave functions in random matrix theory so the statistical theory predicts zero 
correlation between neighboring conductance peaks, but in one of the experiments || correlation was clearly 
present in the form of a slowly varying envelope modulating the peak heights. In subsequent years a number of 
different effects were investigated as candidates to explain this correlation. The simplest is the effect of nonzero 
temperature: since excitation above the Fermi level is possible, several resonances contribute to each peak and 
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a given resonance contributes to several neighboring conductance peaks, inducing correlation. However, in a 
detailed study, this was found to be insufficient to account for the observed correlations JtJ . Other explanations 
that were explored include correlation due to spin-paired levels @||, due to a decrease of the effective level 
spacing found in density functional calculations ||, and due to level anticrossings in interacting many-particle 
systems While these latter explanations rely on subtle electron-electron interaction effects, here we argue 
that peak height correlations already arise within an effective single-particle picture of the electrons in the 
quantum dot. The specific internal dynamics of the dot, even though it is chaotic, modulates the peaks: 
because all systems have short-time dynamical features, chaos is not equivalent to randomness. 

While the statistical theory is "universal" in that it depends on no specific features of the quantum dot at hand, 
the classical dynamics in the dot is clearly not universal. Thus, while correlations between the conductance peak 
heights are generally present in quantum dots, the particular correlations in a given dot are not universal but 
rather involve detailed information about the dot. The simplest information to include is the spatial correlation 
function of the wave functions — this is very short length dynamical information — and an approach including 
this effect was given in Ref. Going beyond this, we use semiclassical techniques to derive a relation between 
the quantum conductance peak height and the classical periodic orbits in the dot. 

The main result is that as a system parameter varies — the magnetic field, for instance, or the number of 
electrons in the dot (controlled by varying a gate voltage) — the interference around each periodic orbit oscillates 
between being destructive and constructive. When the interference is constructive for those periodic orbits 
which come close to the leads used to contact the dot, the wave function is enhanced near the leads, the dot- 
lead coupling is stronger, and so the conductance is larger. Likewise, destructive interference produces a smaller 
conductance. The resulting modulation at frequencies corresponding to the periodic orbits can be substantial. 
Because of dephasing effects, only the short periodic orbits, indeed perhaps only the shortest one, is likely to 
be significant. 

Similar short-time dynamical effects have been noted in other contexts such as atomic and molecular spec- 
tra ]T^|f^] , eigenfunction scarring |l4|,[f5| , magnetotransport in antidot lattices [ fl6| , and tunneling into quantum 
wells [|17|-|2C|] . The periodic orbit modulation that we discuss here is completely omitted in theories in which the 
wave function is assumed to change randomly as the system changes pp|] . Reassuringly, the predicted dynami- 
cal modulation is of the type in the original anomalous experiment 07 More recently, other experimental data 
has been published which show the effect [||,^f) , but to date no systematic experimental study of this effect has 
been performed. 

In the rest of this paper, we generalize some results that have been previously reported in Rcf. j272j to 
address asymmetric lead placement and to incorporate temperature dependence. The derivation given here is 
completely different from the previous one which relied on the methods of Ref. |2(J] : here our approach in terms 
of a statistical ansatz for the wave functions yields more results for chaotic systems but misses the results for 
regular systems that we obtained previously. It has been suggested that asymmetric lead placement would not 
produce an observable oscillation in the average conductance Jl^] , but the method employed there only included 
spatial correlations in the wave functions and not the short-time dynamics which we consider here. In the first 
section we express the height of the conductance peak in terms of the resonant wave function. The basic ansatz 
for the distribution of the wave functions, including dynamical effects, is presented in Section QI. In Section ^ 
our results for the conductance peak heights are obtained. Comparison to numerical results for the stadium 
billiard in Section ^ confirms the adequacy of the semiclassical approach. Finally, we close with a summary 
and discussion of future directions. 



II. THE HEIGHT OF A CONDUCTANCE PEAK IN COULOMB BLOCKADE 

Our starting point is the connection between the Coulomb blockade peak heights and the widths of the levels 
in the quantum dot. This connection is well-known p3[ ; it allows us to express the conductance in terms of 
single-particle quantities. We consider a dot close to two leads so that the width, T, of a level comes from 
tunneling of the electron to either lead. When the mean separation of levels is larger than the temperature T 
which itself is much larger than the mean width, the electrons pass through a single quantized level in the dot, 
and the conductance peak height is |2TJ] 

_ e 2 tt ri r 2 
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where Ti and T2 are the partial decay widths due to the tunneling into a single lead, and spin degrees of freedom 
are neglected. In particular, when the leads are identical and symmetrically attached to the dot, 



Gpcak -XifcT ri (2) 

The partial width is related by Fermi's Golden Rule to the square of the matrix element for tunneling between 
the lead and the dot, M l ~ >d . A convenient expression for the matrix element in terms of the lead and dot wave 
functions, and respectively, was derived by Bardeen [^4j and can be expressed as [ ^7|JT^ 1 

M t-*d = *_ r rfr ^( r)v ^ d ( r ) (3) 

Js 

where the surface S is the edge of the quantum dot. The partial width, T, then, depends on the square of 
the normal derivative of the dot wave function at the edge weighted by the lead wave function. The dot wave 
function ^d in (|J) is calculated for the effective potential, which accounts for interactions in the dot in the 
mean-field approximation. For the partial width we then obtain 

T a [* d ] = ^ V p [ t ] [ dri • V*rf (ri) / dr 2 ■ (r*)* M a) (n)* ^ (r a j\ (4) 
m * e J s Js 1 J 

where a is the index of the lead, the integer £ represents different transverse subbands in the lead, and pe is the 
density of states in the lead for a given subband. To obtain the statistics of the conductance peak heights, we 
thus need to know the statistical properties of the dot wave functions ^>d- 



III. WAVE FUNCTIONS IN THE DOT: THE STATISTICAL DESCRIPTION 

For a single dot, we consider an ensemble of Coulomb blockade peaks — measured either in a narrow interval 
of gate voltage or obtained by following a single resonance under continuously changing magnetic field. The 
wave functions associated with the peaks of the conductance will vary — or "fluctuate" — in a way characterized 
by a distribution P[ip] which we seek. 

It was first conjectured by M. V. Berry that the wave functions of a classically chaotic system fluctuate 
with certain universal properties and can be characterized as random variables |p5| . This is the foundation 
of the first statistical theory of peak heights 0. Subsequently, the statistical ansatz made by Berry has been 
further developed. One direction of refinement is the incorporation of some short length-scale aspects of the 
real classical dynamics. First, a constraint of an arbitrary correlation function 

C{y 1 ,t 2 ) = y^P[V>]<0*( r i)V>(r 2 ) (5) 

was incorporated into the ansatz p6| , p7| . By using the correlation function of a random superposition of plane 
waves, the probability distributions of level-widths and conductance peaks in the case of multi-mode leads to the 
quantum dot were found JH],^] . A distribution similar to this ansatz was derived microscopically for disordered 
systems, a specific kind of chaotic system, using the nonlinear sigma model p8H3lf. 

The next step was to constrain the correlation function by the short-time classical dynamics. Using the short- 
path semiclassical correlation, Srednicki and coworkers | j32|]33f| studied correlations in chaotic eigenfunctions at 
large separations and found that the predicted correlations are in excellent agreement with numerical calculations 
in chaotic billiards (32). This semiclassically constrained ansatz for P[i/j) is much harder to justify — certainly no 
derivation in disordered systems can be made. However, progress towards this goal has been achieved by Kaplan 
and Heller by treating the nonlinear effects of classical recurrences [[l5| . In a recent paper by Kaplan |34| short- 
time dynamics were incorporated into the general probability distribution of Ref. to improve the random 
matrix theory results for the conductance peak height statistics. 

Here we use a maximum entropy technique |35| l to derive the specific form of the distribution P[ip] that 
we need. An advantage of this approach is that arbitrary constraints can be introduced, as in the case of 
normalization which we discuss below. We make the following ansatz: the distribution P maximizes the 
information entropy p& 
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H = - J Dip P [ip] log P [ip] 



(0) 



within the space allowed by the constraints. Here the measure corresponding to the distribution P [tjj] is defined 
in the standard way |33j] 



Vip d = lim U^ =1 dip d (r„) 



(7) 



so that the product P [ipd] Dipd represents the probability that a wave function ip (r) of the original ensemble is 
between ipd (r) and ipd (r) + dipd (r) for any point r inside the dot. 

Assuming that the only constraint imposed on the ensemble of wave functions is the correlation function 
C (ri,r2), the maximum of the functional (Q) under the constraint (|^) is equivalent to the extremum of the 
functional 



F[ip] 



Jv^j -P[ip]logP[iP]- fdvi fdr 2 A(ri,r 2 ){V*(ri)V(r 2 )P[V']-C , (ri,r 2 )} 



(8) 



where the Lagrange multiplier A(r 1; r 2 ) can then be determined from Eq. (||). Setting the first variation of 
F [ip] equal to zero, we find that P[4>] is Gaussian. The final result, obtained by substituting Eq. (||) to find 
A(ri,r 2 ), is 



P [ip] — A exp 



-| J dri J dv 2 r (nJC-^r^raJ^Cra) 



(9) 



where A is the normalization [independent of ip (r)], and G 1 is the functional inverse of the two-point correlation 
function G (r 2 , ri |e) 



Jdr 3 C-^n.rsJCCra.ra) =5(r x -r 2 ) 



(10) 



The coefficient (3 — 1 for a system with time-reversal invariance, when the wave functions can be chosen real, 
and /3 = 2 otherwise. 

It has been shown |?3| that in the small-S. limit for classically chaotic systems, the correlation function 
C (r 2 , ri) can be expressed in terms of the semiclassical approximation to the Green function G sc (r 2 , ri) as 



C(r 2 ,n) 



TPs. 



-ImG sc (r 2 ,ri) + [K 



-3/2 



(11) 



where p sc (e) is the smooth part of the density of states (DOS) in the dot, given by the leading order (Thomas- 
Fermi) semiclassical approximation to the DOS. 

In the semiclassical approximation, the energy-averaged Green function can be expressed in terms of the 
classical trajectories (labeled by the index j) OJ37J, 



G sc (r 2 ,ri) = G (r 2 ,r x ) 



1 1 



in ^/2mh 



^^exp(i| 



inj- j exp 



r?W 2 



2h z 



(12) 



where Sj — Sj (r 2 , ri) is the classical action, Tj is the period, the integer nj is the topological index p2l of the 
trajectory j, and the amplitude Dj is 



Dj = det 



d 2 S j (r 2 ,r 1 ) d 2 S J (r 2 ,r 1 ) 
dr 2 dri de dri 

d 2 S J (r 2 ,r 1 ) d 2 S ] (r 2 ,r 1 ) 
de dr 2 de 2 



(13) 



We have specialized to two spatial dimensions and the last exponential in Eq. (|12| ) is due to a Gaussian averaging 
over an energy window of width W described below. The function Go (r 2 ,ri) is the contribution of the non- 
classical so-called "zero-length" trajectories, those with actions less than or of order h. Therefore, Go (r 2 ,ri) 
cannot be obtained using the stationary-phase approximation, but may be evaluated ]25| , |33]| by replacing the 
actual propagator (r 2 |exp(— iHt/K)] Y\) by its free space analog 
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iHt\ , f dp /.p(r 2 -ri)\ / iH(p,r )t\ 

{ r 2 exp , - _ j ri > « y ^ exp j cxp (-- j (14) 

where ro = (r 2 + ri) /2. The corresponding Green function is then 

n l \ f d P f -P - ( r 2 - r i)\ 1 nK s 

G (r 2 ,ri) = / exp i (15) 



(27rft) 2 V » / e-ff(p,r ) +*0 ■ 

Note that because of the short trajectory involved, this part of the Green function varies very smoothly as a 
function of energy. The smooth part of the correlation function which results is 

Cbte.rO = J- / -^cosf P(r2 ' ri) ) ( 5( £ -g(p ; r )) , (16) 
Psc J (2?m) V n J 

and so Co (r 2 , ri) oc Jo (p |r 2 — ri| /%). This smooth part of the correlation function is rather local in that it 
decays monotonically with separation. Thus, having fully specified the correlation function we wish to use, we 
finally obtain 



P (ipd k) ~ exp 



J dv x J dr 2 r (n)G- 1 (r 2 ,n| £ )^(r 2 ) 



(17) 



A few remarks are required about the width of the energy window W. In the semiclassical limit there arises 
an increasingly broad separation between the short-time dynamics that give rise to system specific behavior 
and the long orbits that are responsible for generating universal statistical fluctuations p!5| . The width W is 
chosen such that the short periodic orbits are included in the sum essentially undamped whereas the long orbits 
are eliminated since their contributions are already accounted for in the statistical ansatz. For the rest of this 
paper, we will eliminate the explicit dependence on W and the sum is understood to contain only the linear 
dynamics. 

The general ensemble defined by the distribution (^) has, however, certain limitations. Strictly speaking, in 
its general form this ensemble is only suitable for calculations of those observables which can be represented in 
terms of only two-point products tp* (ri) ip (r 2 ). The reason for this problem is as follows: instead of the proper 
normalization of each member of the ensemble, 

y dr |V(r)| 2 = 1 , (18) 

the normalization of the wave functions is satisfied only on average, 

Vip(r) P[^(r)] / dr (r)\ 2 =1. (19) 



As a result, the higher order moments, A„ = (J dri . . . J dr n \tp (i"i)| 2 . . . \ip (r n )| 2 )^,, of the distribution are 
different from unity. Therefore, in its general from, the ensemble defined by (^) is not suitable for calculations 
which are sensitive to the n > 1 moments of the distribution P [tp] , such as for the description of the residual 
interactions in quantum dots |38| Eo ]. 

The method developed in this section yields a straightforward way to generalize the distribution (^) to properly 
account for the higher moments. For example, adding an additional constraint 

J V^PUj] J dri J dr 2 |V(ri)| a |^(r a )| 2 =1, (20) 



to the variational problem (|6|) will yield a generalization of the distribution (|17|) which properly accounts for 
the moment A 2 . 

Note in contrast that the errors in the higher moments, n>l, produced by the semiclassical distribution ( |l7j ) 
are of higher order in h, S n ~ O (K 2 ), than the terms taken into account in G sc EjJ. As long as these higher- 
order corrections are not relevant for the quantity under consideration, one can generally use the semiclassical 
distribution (TFi 
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IV. PEAK HEIGHTS DISTRIBUTION 



Since the Coulomb blockade peak heights are uniquely determined by the corresponding dot wave functions 
ipd, the peak heights distribution function P (G) is given by 



P(G) 



Dlpd P (ipd) S(G- Gpeak hM) 



(21) 



where G pca k [ip] is determined by Eqs. (|l|), (||), and (||). The width T depends only on the wave function near 
the boundary of the quantum dot, as follows from Eq. (^). If the function Ps (ip) represents the distribution of 
the wave functions in a narrow strip S along the boundary of the quantum dot, so that 



ip(r), r e S 
j> (r) , r i S ' 



then the conductance distribution is 

P (G) = / Vi> P s [$\e\5(G-G veak 



(22) 



(23) 



The "edge" distribution P$ can be obtained from the general distribution P [ip] by integrating out the values 
of ip, 



Ps [$\ 



Vip P 



tp(tp,i>\ 



As the distribution P [ip] is Gaussian, the resulting functional integral can be calculated exactly, yielding 



Ps 



As exp 



dqi / dq 2 ip* (qi) A^(qi,q 2 )-0(q 2 ) 



where 



#(qi,q 2 ) =G- 1 (q 1 ,q 2 ) 



dq.3 



Q\S 



n\s 



dq 4 G 1 (qi,q 3 )G(q 3 ,q 4 )G 1 (q 4 ,q 2 ) 



(24) 



(25) 



(26) 



and As is the new normalization constant. The spatial integrals are over the part of the total space Q which is 
orthogonal to the edge S, denoted Q\S. 

As follows from Eqs. (pd|)-(12), the "non-diagonal" part of the correlation function is of a higher order in h, 
~ 0(\/h), compared to the "diagonal" part, Go ~ 0(1). The second term in Eq. ( pr^ ) involves the correlation 
functions G (qi,qs) and G (q4,q2), taken between the points of the different parts of the dot, the edge strip 
S for one coordinate and the internal region fl\S for the other. It is therefore of higher order in h, ~ 0(h), 
than the first contribution, G _1 (qi,q2) ~ <5(q2 — qi)O(l) + 0(Vn~). Keeping such higher-order terms is not 
consistent with the leading-order semiclassical approximation we used for G (qi, q 2 ). We therefore obtain 



Ps 



As exp 



2 / dqi I dq2 ^* ( qi ) C 1 ( Cli ' q - 1 ; ' 1 q ' 2 1 



(27) 



An alternative to the argument given here proceeds by noting that integrating out ip should yield a Gaussian 
in ip, and that this Gaussian, by construction of the ensemble, must reproduce the correct two point correlation 
function G(qi,q2). This alternative argument [^7j immediately yields the functional form (p7j). 

When the closed dot is defined by the Dirichlet boundary conditions, the wave function in the narrow strip 
S near the "edge" can be represented as 



ip = zip (y) 



(28) 



where y is the coordinate along the boundary of the dot and z is in the direction of the normal. In this limit, 
the correlation function is 



C(qi,q2) = z 2 z 1 d n C {yi,V2) 



(29) 
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where d n C (2/1,2/2) is defined as the correlation function of the normal derivatives of the wave function at the 
boundary of the dot and can be obtained as 



d n C (y 2 , yi) = -^Im9„G sc (y 2 , yi) + O (h {d+1)/2 ) 



where 



d n G{y 2 ,yi | e) = ^ 



dn^tn (2/2,0) d n i> m (j/l, 0) 



e m - e + iO 

rn 

The semiclassical approximation d n G sc for the normal derivative Green function was derived in Ref. |2C 



(30) 



(31) 



d n G sc (y 2 ,y 1 ) = d n G (y 2 , yi) 
4 



ih 3 \J2-Kih 



(32) 



where fij and [pj] n are, respectively, the Maslov indexes |l2"[ and the normal component of the classical mo- 
mentum of the trajectory j. 

In order to connect the dot wave functions to the lead, let {<p m (y)} be the complete orthogonal set of the 
wave functions corresponding to the transverse potential of the lead. Using this basis, we represent the function 

<p (y) as 



(y) = ^2 a ™§™ (y - vt) 



(33) 



where yi is the contact point of the lead. Assuming that the tunneling between the lead and the dot is dominated 
by the contribution of the lowest transverse subband of the lead and using for the partial width r Q , we 
obtain 



dy </>o (y) (y) 



2nh 



(«) I |5 



(34) 



where is the density of states in the lead corresponding to the lowest transverse subband. For an arbitrary 
moment of the partial width (r™) we therefore find 



(r™> « / da 



2irh (ot) 
mi 



I 1 2m 

| do I exp 



|oo| / dyi / dy 2 <j) (t/i) 



1 



d n C(yi,y 2 ) 



4>o (2/2) 



(35) 



To give explicit expressions for the distribution of level widths and conductance, we specialize to the time- 
reversal symmetric case (/?= 1, GOE) for the rest of this paper; the case when time-reversal symmetry is broken 
by a magnetic field (/3 = 2, GUE) can be treated in an analogous way. In the presence of time-reversal symmetry, 
the wave functions, and hence the coefficient a , can be chosen real, yielding 



(r m ) cx / dr r 



where 



r = 



mi p sc 



dyi / dy 2 (j) (yi) 



t exp(-r/2r) 

Vr 



Im [d n G sc (y 2 ,yi |e)] 



0o (2/2) 



Thus the partial width is characterized by the Porter- Thomas distribution 

1 



P (r) oc —= exp — = 



2r 



(36) 



(37) 



(38) 



with the slowly varying local average T(e). This explicit result for the distribution of level widths is the main 
result of this section. 
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The conductance distribution P (G) can now be simply derived in two limiting cases: (i) when the leads are 
placed symmetrically, so that Ti = I^ [cf. Eq. (||)], and (ii) when one of the partial widths is substantially 
smaller than the other, Ti ^Cl^. In both these cases G ~ Ti [as follows from Eq. (jl])], and the conductance 
distribution is also Porter-Thomas. The "local average" conductance, G, is given by 



G(e) 



2^7™ 



(39) 



where the "local average" width T (e) is defined by Eq. (|37|) , and 7=1 for ri^IV while 7 = 2 for Ti=r2. 

In the general case, Ti/T2 ~ O (1) but not identical, however, an exact calculation of the conductance 
distribution is complicated by the essentially nonlinear dependence of the conductance on the partial widths Ti 
and IV In order to calculate the actual conductance, we choose the area S as the composition of two narrow 
strips, Si and S2, near each of the leads. Using the transverse lead wave functions as the basis in each of the 
two strips, 



(40) 



where the coordinates y^ and yf^ represent the "contact points" of the leads. The partial widths T\ and T2 
are then given by 



2irh (i) 
Fl = —TPq 



r. 



(2) 
2 ^0 



(41) 
(42) 



Assuming equal density of states in the leads, p '(e)~p '(e), for the conductance G we obtain 

G=-— ^pu- 



tt e 2 h 3 i-WnuWia 



2 mikT 



I Op 

.Ml 



-,(2) I 



(43) 



An arbitrary n th moment of the conductance G, (G n ), can now be calculated by integrating over the coeffi- 
cients {a m } for m 7^ 0, yielding 



(G n ) cx / da { Q ] I da 



4 2) 



_(!) 



,(2) 



,(2) 



exp 



„(i) 4 „« 
-a >liia 



«( 2 ) A 

a Q ^22^0 



9„M* A 



(44) 



where the matrix A is 

A a /3 = dyi dy 2 4>o (yi - y| a) ) <9 n G _1 (yi,y 2 ) (po f : 



2/2 



(45) 



Note that the definition ( |45| ) implies that the diagonal elements of the matrix A are proportional to the corre- 
sponding partial widths, An ~ Ti, A22 ~ T2. A straightforward evaluation of the integrals in Eq. ( ff4| ) using 
the substitution 



2 2mlkTG 



Tre 2 h 3 po 



yields 



(G") = y dG G" P (G) 



(46) 



(47) 



where the distribution is 
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P(G) 



Lexp (-Iti AG 
G V 2 



dx 



G 



c 3/2 



■ exp 



1 A 1 A G2 

-~Aux - -A22 — 

2 2 x 



cosh 



A 12 G[J- + - 



x G 
G x 



(48) 



Note that it is only the term involving A12 which makes the remaining integral non-Gaussian and so hard to 
perform. However, this term is semiclassically small: from Eq. (f45| ) it follows that the leading semiclassical 
term in the off-diagonal part of the matrix A is of next order in h compared to the leading diagonal terms. 
The a;-integral in Eq. ( |48| ) is therefore dominated by the interval between 1/ (GAn) and GA22, where the 
off-diagonal matrix element A makes only a small correction quadratic in A\2- Such a correction corresponds, 
however, to higher-order terms in h. Corrections of this order were already neglected in the original semiclassical 
expansion of the Green function, and so to be consistent we discard all effects of the off-diagonal matrix element 
A12 here. The integral in Eq. ( [f8|) can now be easily performed. 

The semiclassical approximation to the conductance distribution is, then, simply a Porter-Thomas distribu- 
tion, even in the general case: 



P{G) = 



1/2 



1 



exp 



G_ 



(49) 



where the "local average" conductance G is 




Gee 



and the "partial conductance" G a is related to the partial width r Q via the standard relation 

G a ee (Tre 2 /2hkT)T a . 



(50) 



(51) 



As the semiclassical Green function, G sc , and, consequently, the correlator, d n C, can be expressed as a sum of 
the contributions of "zero-length" and longer classical trajectories, similar decompositions hold for the average 
partial width and average partial conductance: 



— osc 

*• rv 



r — r 

— — — o 

G a = G a + G p 



(52) 
(53) 



where the "oscillatory" parts, V and G , depend on the longer classical trajectories and are of next order 
in h compared to the smooth contributions, T and G which are the zero-length contributions. A consistent 
semiclassical approximation, as in Eq. (|l2]), then requires expanding G and keeping only the linear terms in 
the oscillatory contribution. We thus obtain 



G 



5! 



G9 



t 



gT 



1 



1 + JcfM Gi 1 + Jc?M 



Gn 



(54) 



We now proceed to the semiclassical calculation of the "local average" partial width T. The defining equation 
( [37| ) involves the functional inverse of the Green function, which is a hard object to calculate. Instead, we will 
use the original definition (Q), which for the local average partial width yields 



- 27r7T 



dyi / dy 2 cj> a (yi - ye.) <t> a (j/2 - Hi)* d n C (yi,y 2 ) 



(55) 



where the correlation function of the normal derivatives of the dot wave functions d n C (yi, 2/2) is related to the 
semiclassical Green function by Eq. (p0[). 

If we now use some information about the lead wave functions, we can obtain an explicit expression for the 
average width T in terms of the classical dynamics in the dot. When, as we assumed above, the tunneling from 
the lead to the dot is dominated by the lowest transverse energy subband in the constriction between the lead 
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and the dot ||, the transverse potential in the tunneling region can be taken quadratic: Ue ~ K (y — ye) ■ In 
this case, the transverse dependence of the lead wave function is simply a harmonic oscillator wave function, 
so that at the edge of the dot (j>o — ce exp[— (y — ye) 2 /2a 2 s ], where ye is the center of the lead and constriction, 
and the effective width is a c fr = Vfr/ \/2nm*. While the exact form of the lead wave function is not crucial, the 
^-dependence of the width is important for the semiclassical argument which follows; note that a e ff ~ \fh. does 
not depend on a particular transverse potential. 

Using this information about 0o in the expression for the diagonal matrix elements An, A22, we see that the 
lead wave function restricts the integration to a semiclassically narrow region of width a e e ~ y/h. This allows 
one to express the contribution of the open trajectories entering the Green function in terms of an expansion 
near their closed neighbors: 



r = r o + / dy 

m 



dpyfw (y,Py)^2 



x exp 



.S a (y, 0; y, 0; e) 



to?! + mf 2 + 2 



exp 



1 

h to? 



2m? 2 



(56) 



where Tq is the monotonic part of the resonance width, (pi) n and {pf) n are the normal components of the initial 
and the final momenta of the closed orbit a, the momentum p = (p* + p/) /2, the 2x2 monodromy matrix 
p2| M a = (mf ■ ) is defined via the linearization of the Poincare map near the closed orbit a , and calculated 



at the contact point near the lead. We have also introduced in Eq. 
wave function 



(|56|) the Wigner transform of the lead 



f^iy^^h- 1 / dAy 6, (y- Ay/2,0) </>S (y + Ay/2, 0) exp (i Py Ay/h) 



(57) 



which describes the distribution in transverse position and momentum of electrons tunneling into the dot. 

In leading order in the distance between the contact point y of the closed orbit a, and the center of the lead 
yi, the action of the closed orbit S a scales linearly: 



S a (y, 0; y, 0) cx S ( W , 0; y t , 0) + Ap« (y - ye) 



(58) 



where Ap y is the change of transverse momentum after the traversal of the closed orbit. Assuming e.g. a 
Gaussian form of the lead wave function, the contribution of each of these closed orbits is suppressed by a factor 
exponentially small in Ap y . This suppression is the effect of the mismatch of the closed orbit (momentum) with 
the distribution of transverse momentum at the lead, which is centered at zero with width 5pe ~ h/a e s ~ \/% for 
the lowest subband. Therefore, only closed orbits with semiclassically small momentum change Ap contribute 
to the width. This in turn implies that the closed orbit is located semiclassically close (within a distance ~ Vft) 
to a periodic orbit for which Ap = 0. Using this proximity to a periodic orbit we can re-express the actions and 
momenta of the injection orbits in terms of the properties of their periodic neighbors (labeled by the index fi) 
as follows: 



S a (y» 0; y, 0) ~ Sf,, + 



Tr [M p 



2^12 



— {y~y^) 2 i 



'11 



'22 



(y - Vn) 



(59) 
(60) 



Substitution of (&9) and (pfl) into Eq. (|5q) and integration over y yields 



r = r + ^ cos if 



fi-.p.o. 



Si, 



(61) 



where the monotonic part is 



r = ^c|a eff ^ e-C [7 (C)+/i(C)]. C = 

TO* 



2h 2 ' 



(62) 



the amplitude is 
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h rfPz r T 2,, f W1 , _2W, , J\rl/1 



A, = 4^ ^ [Tr 2 [M„] (1 
(1+4) 




x exp( - , a + p2 9 , - "7%, ) (<)>!) 



ith 



TOi2-m 2 i ± 



^/ (m 22 -mii) 2 + (m 2 i+mi 2 ) z 



2m y / ", (i 



J Tr [M M ] + 2 V ?i 

fl^arctanfg^g^ (64) 

y = cos6 (y M - y^) /a off + sin0 p^a cff /h 
p = cos 8 p%a e g/h - sin0(y M - yt) /a c g , 

and (f>fj, is a slowly varying phase. Here /„ is the Bessel function of complex argument, p is the magnitude of the 
electron momentum, is the electron momentum for the periodic orbit /x at the bounce point (turning point), 
y^ is the bounce point coordinate, is the action of the periodic orbit, and = (rnfj) is the corresponding 
monodromy matrix JlJ] . Note the sharp suppression of the oscillatory effects in Eq. ( |63|) if the periodic orbit does 
not match up to the lead wave function in both position and momentum space. The mismatch is characterized 
by y and p; the most favorable case is that of a perpendicular periodic orbit hitting the edge of the dot right at 
the center of the lead, p^ = and y^ = yt so that y=p = 0. 

An explicit expression for the average conductance follows from the relation between the partial width and 
the partial conductance Eq. j5l|). Using Eq. ((53), we see that G can be written in the form 



G = Go + B » cos ( -T + <t>» ) ( 65 ) 



\~" n I A* 

=0 , ^=;0 



where is simply related to A^, G 1; and G 2 - This together with Eq. ( |49| ) defines both the average conductance 
and its fluctuations. 

Further characterization of the peak fluctuations can be obtained from the peak-to-peak correlation function: 
this is a particularly interesting quantity because of the correlations sometimes observed experimentally [||J|] , 
as discussed in the Introduction. A natural measure of the statistics of nearby peaks is given by 8G (E m ) = 
G (E m ) — (G (E n )) n in terms of which the correlation function is 

Corr m [SG, 5G] = (SG (E n+m ) SG (E n ))J ([6G(E n )] 2 ) n . (66) 

Substituting the conductance distribution (^) into (|66|), we obtain 

4G 2 + 3E M ^ 



Corr m = S m>0 + (1 - 8 m>0 ) x — ,^ 2 , ov ^ n9 — ■ (67) 



Throughout this paper we have concentrated on energy (or equivalently peak number) as the tuning parameter 
causing the peak height variation. This is just an example: exactly analogous considerations apply to any 
parameter causing changes in the wave functions of the quantum dot. In particular, similar oscillatory behavior 
is expected in the height of a given peak as a function of magnetic field, often the most experimentally accessible 
parameter. As the field varies, the change in the action of a periodic orbit is proportional to the (directed) area 
that it encloses. Thus, the peak heights should exhibit an oscillatory envelope whose frequencies are proportional 
to the areas of the periodic orbits. 

V. COMPARISON WITH NUMERICS AND EXPERIMENT 

Since one of the main theoretical results of the present paper concerns the periodic modulation of the Coulomb 
blockade peak heights, it is natural to consider the Fourier power spectrum of G pea k(fc). In Fig. |l| we present a 
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FIG. 1. The peak conductances (left column) from tunneling through subsequent energy levels in the stadium 
quantum dot and the corresponding "length spectra" / (L) (right column) for different lead configurations (shown in 
the insets). In the peak conductance plots, each peak is placed at the wavevector k corresponding to its level; R is the 
radius of the half-circle parts of the stadium dot. A Gaussian lead wave function appropriate for tunneling from a single 
transverse mode is used with width ka c a — 15. The red curves represent the semiclassical envelopes, defined by the 
contributions of the relevant periodic orbits (top: "diameter" orbit, middle: "V'-shaped orbit, bottom: both diameter- 
and V-orbits). Length spectrum of the oscillations in G(k) obtained from the Fourier power, numerical (green) and 
semiclassical (red) results are compared. The power is normalized to the mean conductance. The arrows at the top show 
the positions of the relevant periodic orbits and their repetitions (black arrows), as well as the "combination lengths" 
L\ ± L2 (blue arrows). In the top panel, the peak at L/R ~ 4 is the diameter, that at 8 is its repetition. In the middle 
panel, the peak at L/R = 2 (l + \/2) ~ 4.8 corresponds to the V-shaped orbit, the peak at L/R ~ 9 represents its 
repetition. In the bottom panel (asymmetric leads), the broad peak at L/R ~ 4.5 represents the total contribution of 
both diameter- and V-shaped orbits. For the stadium dot, the principal peaks appear at 4 and 4.8, because we use only 
the wave functions symmetric about the vertical and horizontal symmetry axes (equivalent to using only the even-even 
states of the stadium). Note the excellent agreement between the semiclassical theory and the numerical results for 
symmetric leads, and adequate representation of the principal peak for asymmetric leads. 

comparison of the numerical and semiclassical power spectra, calculated for a chaotic (stadium) dot, for three 
different placements of the leads. Previously, we reported the case for leads placed symmetrically as in the 
upper plot of Fig. [I] for kR — 70. The exact conductance peaks are obtained numerically from Eq. (|l]) with 
the eigenstates being constructed using the the method of Ref. [Q. To observe the variation in peak height, 
we vary the energy, or equivalently the wavevector k — p/h, which changes the number of electrons on the 
dot as more levels are filled. The data clearly demonstrate that the power spectrum has well-defined peaks 
corresponding to periodic orbits. The numerical results for the symmetric leads show excellent agreement with 
the semiclassical prediction. 

The situation is however different for asymmetrically positioned leads when there is no single short periodic 
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FIG. 2. The dependence of amplitude of the length spectrum peak at L/R ~ 4 on ka e fj for kR = 70. The leads are 
symmetrically attached to the middle of the semicircle segments of the stadium dot. 

orbit connecting both leads. In this case, only the main peak corresponding to the first repetition of the 
relevant periodic orbits, the "diameter" and the T^-shaped orbit, is adequately reproduced. The higher-frequency 
behavior, however, is substantially different from the semiclassical prediction. We attribute this difference to 
the non- linear mixing of the oscillations of different partial widths, neglected in our derivation of Eq. (|54|). The 
pronounced peak at the difference length Ly — Lp, where Ly and correspondingly represent the lengths of 
the 1^-shaped and diameter orbits, strongly indicates that, although semiclassically small, the mixing effects of 
higher order terms in Eq. ( |50|) can be significant in the experimentally relevant parameter range. We numerically 
verified that the sum and difference lengths can be partially obtained by Eq. (|50|). 

As follows from Eq. (|6l|), the oscillatory component of the "local average" conductance and the height of 
the corresponding peak in the power spectrum depends nontrivially on the position and the width of the lead. 
This dependence is illustrated in Fig. E| where we plot the amplitude of the "diameter" orbit contribution 
to the conductance as a function of fca ff extracted from numerical length spectrum and the corresponding 
semiclassical prediction. 

In Fig. H| we compare the semiclassical correlation function with numerical data for the stadium dot. The 
oscillatory behavior for large separations reflects the peak in the corresponding power spectrum in Fig. ^ and 
is in agreement with the semiclassical result. The positive correlation for nearest neighbors is also in agreement 
with the semiclassical theory, demonstrating the influence of dynamics even in this apparently non-semiclassical 
regime. 

When T> A, the major source of correlations between neighboring peaks is the joint contribution of several 
resonances to the same conductance peak 0. In this regime the "nearest-neighbor" correlator is Corr m=1 ^l, 
and the dynamical effect accounts for only a small correction to the correlation function. However, for low tem- 
perature T< A, the correlations due to temperature are exponentially suppressed. In this regime, as illustrated 
in Fig. |J, the correlations induced by dynamical modulation dominate, and they account for the experimentally 
observed enhancement of correlations at low temperatures |^] . For finite temperature each resonance is weighted 
by combinations of Fermi-Dirac functions and occupation numbers |23]| . The occupation numbers used in Fig. 
^ were obtained by employing a recursion relation [^3| ; see Appendix |a]. As the temperature increases more 
resonances contribute to a single conductance peak, and thus dampening the effects of the longer orbits. 

In Fig. |E| we present the results of the calculation of the probability distribution of G pca k for a stadium 
quantum dot for both the "symmetric" and "asymmetric" placement of the leads. For comparison, we show 
both the actual distribution, Eqs. ( [49|) (|54|), and the standard Porter-Thomas result without any account of the 
modulation of the average conductance: -P(G poa k) = -\/4/7rG P e a k exp(— G pca k). As the individual peak-height 
distribution is essentially a local measure, it is not strongly sensitive to the correlations, and both the standard 
and the dynamical theories predict nearly the same result, and both are consistent with numerical calculation. 
This explains why no dynamical effect was observed in the experimental peak-height probability distribution 

ill- 

In contrast, the periodic modulation of the peak heights has been observed in several recent experiments 
P,ppl|]. The clearest observation is in Ref. ||: the data in their Fig. [I] show modulated peak heights as a 
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FIG. 3. The peak-to-peak conductance correlation function for (a) symmetrically placed leads (attached to the "di- 
ameter" of the stadium dot), and (b) asymmetric leads (as in the inset to the length spectrum at the bottom of Fig. |lj. 
The numerical correlation function (circles with typical error bars) — the average of all pairs of peaks m peaks apart — is 
in good agreement with the semiclassical theory (red). The agreement for small m is surprising since this regime is 
not semiclassical, but shows how dynamics can give rise to correlations even between nearest-neighbors. The difference 
between the periods of the modulation in (a) and (b) is accounted for by the difference in the values of kR used for 
the correlation function: the calculation for (a) is performed near kR = 70, while (b) corresponds to an interval near 
kR = 140. 




10 20 30 40 

in 

FIG. 4. The peak amplitude correlation function for stadium-shaped quantum dot with symmetrically attached leads, 
for the temperature T = 0.25A(red curve), 0.5A(blue), A(green), 2A(yellow) and kR = 70. 

function of the number of electrons in the dot. In their trace of 90 peaks, approximately six oscillations are 
visible, yielding a period of ~ 15 peaks. In our treatment, this period is related to the period of fundamental 
oscillation in Eq. (|6l]) . A variation in action AS^ can arise from changes in either or both the Fermi momentum 
and the dynamics, i.e. the lengths of the orbits. If only the momentum varies, then the fundamental period 

is given by (k^§f-J = h/ T ti where is the period of the relevant orbit, and the ratio of this to the level 
spacing A gives the period of the peak heights. In the billiard approximation, = L^/vp, where L M is the 
length of the periodic orbit and Vf is the Fermi velocity, which can be calculated from the experimental density 
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FIG. 5. Conductance statistics: probability distribution function for (a) symmetric leads at kR = 140 and (b) 
asymmetrically placed leads at kR — 70. The numerical probability distribution (histogram) is for the entire range of 
data in Fig. ^ and is compared to both the semiclassical theory (red) and the standard statistical theory based on random 
wave functions (green) . The two theories predict nearly the same result for this quantity (especially for asymmetric leads, 
where the dynamical modulation is weaker), and both are consistent with the numerics. 



p4j . Using the appropriate spin-resolved level spacing A ~ 10 fj,eV (which is half of the spin-full value from 
the measurements in Ref. j8j), the experimental period implies L„ ~ 4.5/ito. This value is inconsistent with the 
period given by the shortest periodic orbit contacting the leads |4f| — the orbit from the lead to the "pin" gate 
and back whose length we estimate to be 0.9/im. 

Changes in the dynamics should also be considered. By examining the configuration of the dot in the insert of 
Fig. 1 of Ref. ||] , the gate voltage appears to be situated on the shortest periodic orbit of each lead. By making 
the gate voltage more negative the electron will have a shorter path and thus contribute to the change of the 
action in Eq. (]6l]). If the Fermi energy of the system remains constant, then we can calculate a plausible range 
for the period using two extreme simplified models for the deformation of the boundary. First, the gate voltage 
is modeled as a small local semicircular deformation. Equating the area of the semicircle to the number of peaks 
times the change of area caused by adding one electron on the dot without a change in the Fermi energy, we 
estimate the period to be ~ 3 peaks independent of the Fermi energy. The other extreme is considering the 
entire side to move uniformly. The same procedure yields a period proportional to the square root of the number 
of electrons on the dot. Assuming that the typical dot in these experiments has 100 electrons ||, we obtain a 
period of ~ 13 peaks. Thus, the experimental result of a period of 15 peaks cannot be obtained solely by the 
orbit length change, but must also include a change due to the difference in Fermi energies of successive peaks. 
A more detailed model of the gate voltage is necessary in order to make a better prediction of the oscillation 
period. 

A similar approach to the peak modulation as a function of magnetic field is also in agreement with the 
experimental results |||lj , where a quasi-periodic modulation of the peak heights was observed with the period 
AB ~ 35mT. In our treatment, this period is given by the ratio of flux quantum hc/e to the area Ao enclosed 
by the periodic orbit. From the experimental oscillation, we obtain Aq ~ 0.12/im 2 . This is consistent with 
the total area of the dot, 0.32^m 2 ||, considering that there is likely to be some cancellation of fluxes between 
different parts of the orbit. 

A puzzling feature of the initial experiments was that the dynamical modulation of the Coulomb Blockade 
peak heights was not seen in the experiment of Ref. j5j . We attribute this behavior to two factors: the positioning 
of the leads relative to the gate and a relatively small mean free path. First, if the gate used to change the 
number of electrons is not along the shortest periodic orbit of either lead and the Fermi energy does not change 
appreciably in the dot, then one should not observe oscillations in the conductance peaks. In the geometry of 
Ref. H the leads and gate seem to be rather disconnected so this may be a factor. Second, in this experiment 
the mean free path I ~ 0.4/im only marginally exceeds the typical size of the dot d ~ 0.25/xm, while the length 
of the shortest periodic orbit is at least twice the effective "diameter" d of the dot: L min > 2d ~ 0.5/im > I. 
If the mean free path is caused by short-range diffractive scattering, the dynamical effects are suppressed and 
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will not affect the Coulomb Blockade measurements. However, in the opposite limit of a smooth scattering 



potential, dynamical effects caused by coherent branched flow 46 may still be present. The short mean free 
path measured in Ref. || suggests the presence of impurities in the two-dimensional electron gas layer leading 
to a short-range scattering potential and so suppression of dynamical effects. 

The agreement of our semiclassical theory with experiment may seem surprising, since the adding of elec- 
trons changes the effective potential defining the dot because of the added charge. However, experiments on 
"magnetofingerprints" of the peaks p7|] suggest robustness of the effective potential — its change from peak to 
peak seems to be small in this case. In contrast, to affect the dynamical modulation one must substantially 
change the action of the shortest periodic orbit, which typically requires a much larger change in potential such 
as could be caused by the external gates. 

VI. SUMMARY 

In conclusion, using semiclassical methods, we developed a dynamical statistical theory of Coulomb blockade 
peak heights in chaotic quantum dots. We derived the peak height distributions and the correlation functions, 
and showed that the corrections to the corresponding results of the standard statistical theory can be expressed 
in terms of the classical periodic orbits of the dot. Both our analytical results and numerical simulations clearly 
demonstrate that the dynamical effect is significant for both symmetric and asymmetric lead placements. 

We close with two further experiments suggested by our results. First, if the tuning parameter used to change 
the number of electrons, such as a gate voltage, does not change the action of the dominant periodic orbit, then 
no modulation connected to that orbit should be seen. In particular, gates which affect different parts of the 
dot may produce different oscillatory behavior. Second, several samples made in a robust geometry — a circle 
with directly opposite leads, for example — should show the same modulation. Any deviations from the same 
behavior would be a sensitive indication of the material quality. 
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APPENDIX A: TEMPERATURE CALCULATIONS 

For nonzero temperatures the conductance is obtained from a weighted sum over the zero temperature partial 
widths T\ P3| . For symmetric leads this yields 

o-TVrZ-** (A1) 

If kT, A <C e 2 /C, then the weights are given by 

w x = Af(AF No - E F )(n x ) No [l - f(E\(N) - E F j\ (A2) 

where AFn is the change in the canonical free energy from N — 1 to N, (nx)N is the canonical occupation, 
Ef = En + (N— l/2)e 2 /C is an effective Fermi energy and /(e) = [l + cxp(e/fcT)]~ 1 is the Fermi-Dirac function. 

To obtain the canonical free energy and canonical occupation number we use a recurrence relation developed 
by Brack, Genzken and Hansen |^3| for the partition function Z(N,M; /?); N is the number of particles, M is 
the number of levels and (3 = 1/kT. The final result for the partition function will not numerically depend upon 
M for large M. The partition function is formally given by 

Inm 

Z(N, M- f3)=J2 exp(-/?£ Q (A0) = exp(~PE )z(N, M; (3) (A3) 

Q = l 



1G 



where 



In m 



Ct = l 



(A4) 



Here E a (N) is the sum of the energy of the single particle occupied levels em which does not include the 
charging energy, Eq will be defined below, and Inm is the number of ways to fill M levels with N identical 
particles. The recurrence relation derived in Ref. |43] is 



Z(N, M; (3) = Z(N, M - 1; j3) + exp(-/3e M )Z(7V - 1, M - 1; (3) for N > 1, M > N 
with the constraints 

Z(0, M; (3) = 1 VM > 
Z(N,N- l;(3) = VN > 1 . 

Note that the same recurrence relation also holds for z(N, M; (3). The choice E {N) — J2m=i Cm yields the 
result 



(A5) 



(A6) 
(A7) 



z(N,N;/3) = l 



(A8) 



Using the conditions (A7) and (A8) as a starting point for the recurrence relation Eq. ( A5), we obtain 
z(N, M — > oo; (3) and thus Z(N, M — > oo; (3). For the small temperatures that we consider, the convergence of 
the recurrence relation is rapid. 

Similarly, one can calculate a modified partition function Z[(N,M; 0) which has level A removed from the 
spectrum. The probability for level A to be unoccupied, P{n\ = 0}, is, then, simply Z'^/Z. In terms of this 
probability, the average occupation numbers are given by (tix}n = 1 ~ P{n\ — 0}. Finally, the canonical free 
energy for N electrons, F(N), appearing in Eq. (A2) is 



F(N) = - - In Z(N, M -» c5o; (3) 

P 



(A9) 
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